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Abstract 

We examine the detectability through indirect means of weakly interacting dark 
matter candidates that may constitute not all but only a subdominant component 
of galactic cold dark matter. We show that the possibility of indirect detection of 
neutralinos from their annihilations in the Earth and Sun is not severely hampered 
by decreasing neutralino relic density. Upward-going muon fluxes in underground 
detectors from neutralino annihilations in the Sun can remain above the threshold 
of detectability of 10 muons/km 2 /yr for neutralinos composing 1% or more of the 
halo dark matter. Similarly, signals from neutralino annihilations in the Earth can 
also remain high for neutralino densities of 1% of the halo and actually would only 
be observable close to this low density for neutralinos lighter than 150 GeV. We also 
show that there are many models which simultaneously have high direct and indirect 
detection rates making some model discrimination possible if a signal is seen in any of 
the current dark matter searches. 
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We examine the observability of collisionlcss cold dark matter (CCDM) by 
indirect detection when it is merely a subdominant component of the cold dark 
matter. Originally wc were led to consider such a question by claims of the 
need for collisional cold dark matter [1] as the main form of dark matter in 
the universe, but the question is valid independently of such claims. Even if 
non-CCDM is proven to be unnecessary, there remains the possibility of the 
CDM consisting of several populations, the one we are searching for not being 
the dominant one. So the question is: if the previously favored CDM candi- 
dates, such as Weakly Interacting Massive Particles (WIMPs), constitute only 
a fraction, say 1% or less, of the local dark matter halo, would these particles 
still be observable in the current and proposed direct and indirect dark matter 
searches? We could even reverse our question in the following manner: if we 
see a CDM signal in any of our searches, could we be observing a subdominant 
component of the total CDM? Additionally, could a combination of direct and 
indirect detection data resolve the issue of which component of the total CDM 
we are observing? 

Unless there is segregation for different types of cold dark matter, the ratio 
of collisionless CDM to total CDM should be the same locally in the Galaxy 
and globally in the whole Universe. This is expected if all CDM components 
are collisionless. Thus in the following we assume that the local fraction of the 
particular collisionlcss CDM Jccdm is related to the CCDM density parameter 
Qccdm (the CCDM relic density in units of the critical density) through 



Here pccdm is the local density of a particular CCDM candidate, such as a 
particular WIMP, p\ oca \ ~ 0.3 GeV/cm 3 is the local halo density (at the location 
of the Earth), flee dm is the relic density of our particular CCDM candidate, 
and CIcdm — 0.3 is the total contribution of DM to the total energy density of 
the Universe. 

Naively, one may expect that if the local CCDM density is, say, fee dm = 1% 
of the local halo density, the fluxes of particles produced by WIMP- WIMP 
annihilations in the halo and elsewhere, being proportional to some power n of 
the local WIMP density, would decrease by an amount fg c DM , i.e. by 0.01" in 
our example. However, the relic density of thermal WIMPs Q x is approximately 
determined by their annihilation cross-section a a through the relation 



where (a a v) E v is the thermal average of the annihilation cross section times the 
relative velocity of the WIMPs at freeze-out, in the early universe, and h is the 
reduced Hubble constant, h ~ 0.7. Hence a reduction in the relic WIMP density 
requires an increase in their annihilation cross-section in the early Universe. In 



fcCDM = 



Pccdm _ Qccdm 

Plocal QCDM 



(1) 



n x h 2 



3 x 10 27 cm 3 /sec 



(2) 
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general, (<j a v) E v is related to (<r a v), the annihilation cross section of WIMPs in 
the galactic halo or at the center of the Sun or of the Earth. These two averaged 
annihilation cross sections differ in the temperature of the WIMPs when they 
annihilate: about a twentieth of the WIMP mass in the early universe, practi- 
cally zero in the halo, the Sun and the Earth. However, there are exceptions 
to this relation: when coannihilations [2] are important the WIMP population 
has to be followed together with the population of other particles with which 
it can be interchanged. As a consequence, (a a v) E v is an effective annihilation 
quantity that includes annihilation of and with other particles, and which is not 
directly related to (a a v) at zero temperature. Since the annihilation signals are 
in general proportional to the product of the annihilation cross section and the 
square of the WIMP density, and the latter decreases linearly with the annihi- 
lation cross section, a more educated statement is that the annihilation signals 
decrease linearly with WIMP density. This is the case if (<J a v) E v and (cr a v) 
are proportional. We will see below (Fig. 3) that this linear decrease is true for 
the envelope of the largest signals of neutralino annihilation in the halo. 

The situation is more complicated when the signals depend not only on 
the annihilation cross section a a , but also on the scattering cross section a s of 
WIMPs with nucleons. Besides direct detection, this is important for WIMP 
annihilation in the Sun and the Earth. The number of WIMPs accumulated in 
those bodies is proportional to (J s pccDM, i-c. to ascr^ 1 . So at short times the 
annihilation signal goes as a a (a s (T~ 1 )' 2 . At times much longer than the equili- 
bration time between WIMP capture and annihilation, instead, the annihilation 
signal depends entirely on the product a s pccDM ~ ffs^a 1 , since in equilibrium 
the annihilation rate equals the capture rate. 

Now, the scattering cross section is not independent of the annihilation cross 
section. As argued in earlier papers [3, 4, 5], when the annihilation cross sec- 
tion grows, the scattering cross section may grow by the same factor. If so, 
the product <j s pccdm ~ ffs^ 1 remains approximately constant, and signals 
proportional to this product, such as those from equilibrium annihilation in the 
Sun and Earth, remain nearly unchanged. 

Scattering and annihilation cross sections of WIMPs are in general related 
by a crossing symmetry. Unpolarized cross sections are spin sums of amplitudes 
squared, and a crossing symmetry relates the amplitudes of two processes that 
differ by having a particle in the final state exchanged with a particle in the 
initial state, both becoming their respective antiparticlc. For example, the am- 
plitude for neutralino-quark scattering xi — > XI i s related to the amplitude for 
neutralino-neutralino annihilation into quark-antiquark pairs XX ~~ > <M- Often, 
other annihilation channels besides qq are open, such as annihilation into W 
or Z boson pairs, and since these particles are not present in the nucleon, the 
crossing symmetry is imperfect. Moreover, the scattering and annihilation cross 
sections are computed at different kinematical points, and may differ in the spin 
structure, the presence of resonances, the integration limits over momenta, etc. 

The relation between the WIMP annihilation cross section and the cross 
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section for WIMP scattering off nucleons is further complicated by the fact 
that annihilation can be considered to produce quark- antiquark pairs, while 
scattering occurs with quarks bound in nucleons. The annihilation amplitudes 
are functions of the quark masses, while the scattering amplitudes are functions 
of the nucleon mass. Moreover, the scattering amplitude contains the matrix 
elements of the quark currents in the nucleon, which are of course not present 
in the annihilation amplitude into quark- antiquark pairs. 

Despite these imperfections in the crossing symmetry, we can see it at work in 
the case where WIMPs are neutralinos in the Minimal Supersymmetric Standard 
Model (MSSM), with some sets of parameters. We use the version of the MSSM 
described in [6], whose seven input parameters are given at the weak scale, 
as implemented in the DarkSUSY code [7]. In Fig. 1 we make the crossing 
symmetry explicit in a particular example. 

We present a set of models obtained by modifying only the neutralino compo- 
sition (and thus its couplings) starting from a single original model that exhibits 
high muon fluxes from neutralino annihilations in the Earth at low neutralino 
relic density (we argue later that this high flux is due directly to the compen- 
sation provided by the crossing symmetry). The particular parameter values 
of the original model (one of the points in Fig. 1) are as follows: the mass of 
the neutralino m x = 147.9 GeV, the mass of the CP-odd Higgs boson itia = 
113.1 GeV, the ratio of Higgs vacuum expectation values tan/3 = 47.52, the 
sfermion mass parameter too = 1873 GeV, and the trilinear couplings in the top 
and bottom sectors A t = -1.758too and A], = -2.152to . This particular original 
model has a neutralino relic density of Vl x h 2 = 2.9 x 10~ 3 , a neutralino gaug- 
ino fraction Z g — 0.2511, and muon fluxes from neutralino annihilations in the 
Sun and Earth of <E>® = 143.3 muons/km 2 /yr and <I>® = 44.25 muons/km 2 /yr, 
respectively. In order to examine how the scattering cross section changes with 
the annihilation cross section, we varied the higgsino and gaugino mass pa- 
rameters of the original model, /i and Mi, while keeping the neutralino mass 
and all the other MSSM parameters fixed. This variation produces a "line" in 
the model parameter space with constant neutralino mass but varying \i and 
Mi parameters, and hence different gaugino and higgsino compositions. Fig. 1 
shows plots of the spin-independent neutralino-proton scattering cross section 
er x _p, the neutralino- neutralino annihilation cross-section at zero temperature 
and momentum (a a v), and the neutralino relic density fl x h 2 , all plotted versus 
Zg/Zh = Z g /{1 — Zg), the ratio of gaugino and higgsino fractions (defined such 
that Z g + Zh = 1). The range of /i and Mi used to produce Fig. 1 was 150 GeV 
< (j, < 1.2 TeV and 8.9 TeV > M 2 > 290 GeV. 

From the plot we see that for mixed and gaugino-like neutralinos (Z g > 
O.lZh) there is an excellent correspondence between fl x h 2 and (a a v): as the an- 
nihilation cross section in the Earth or Sun increases, the relic density decreases 
and vice-versa, as expected. Furthermore, we see evidence of the crossing sym- 
metry at work by examining the scattering and annihilation cross sections: as 
the annihilation cross section increases, so does the scattering cross section. 
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This correspondence makes neutralino models with low relic densities but large 
direct and indirect detection rates possible. Fig. 1 shows that the crossing 
symmetry is not effective when the neutralino has a small gaugino component 
(Z g < O.lZ/j). One reason the crossing symmetry is not perfect at low Z g is be- 
cause annihilation into W and Z pairs dominates for these neutralinos. Hence, 
the annihilation cross section remains large, but the scattering cross section de- 
creases at low Z g . Moreover, at low Z g , the neutralino is almost a pure higgsino 
and is almost degenerate in mass with the second lightest neutralino and the 
lightest chargino. Therefore, coannihilations between neutralinos and charginos 
make {<J a v) E v larger than (a a v). In fact, taking coannihilations out, we find 
that the neutralino density without coannihilations, shown in Fig. 1 with the 
line labelled "w/o coanns", would be inversely proportional to (a a v). 

In the rest of this paper, we study the concrete case of the lightest neu- 
tralino. For this purpose, we use a table of models allowed by all accelerator 
limits (but without imposing any constraint from the recent measurement of 
the muon anomalous magnetic moment [8]). The table has been produced with 
the DarkSUSY code [7] over the last few years for other purposes, i.e. having 
in mind other issues which were addressed in the papers of Ref. [6, 9] for which 
the models were originally computed. Thus, to produce all the remaining fig- 
ures, we have not done any particular sampling of the models to favor lower 
densities and higher detection rates. We restricted our attention to models with 
tt x < Udm = 0.3 (Vl x h 2 < 0.15) for which we found about 45,000 points in 
parameter space. 

We use this table of models to show the validity of Eq. (2) . In Fig. 2 we plot 
il x h 2 versus the zero momentum (and zero temperature) limit of the neutralino- 
neutralino annihilation cross section times relative velocity (<r a v), together with 
the approximate relation in Eq. (2) (straight line). From the plot we see that Eq. 
(2) is a good approximation for the models with the largest (a a v). Generally, 
points are to the right and below the line because the cross sections increase 
with energy and more annihilation channels, including coannihilations, are open 
in the early universe. Some points are to the left of the line due to resonances 
active at low energy, or above the line due to coannihilations with particles 
which have a lower annihilation cross section than neutralinos. 

In a previous paper [4], the same table of models was used to study the 
possibility of direct detection of a subdominant component of neutralinos and it 
was concluded that many models of subdominant neutralinos, even with WIMP 
halo fractions as low as 10 -4 , are within the discovery limit of proposed detectors 
in the conceivable future. 

We now examine indirect detection of a subdominant neutralino component. 

A proposed method of indirect detection consists of searching for gamma- 
rays and rare cosmic rays produced by WIMP annihilations in dark halos [10] 
or at the galactic center [11]. The rate of these annihilations scales as 

T ~ o aP \. (3) 
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Since, as discussed above, p x scales as fl x and a a scales approximately as l/r2 x , 
the annihilation rate in the halo scales approximately as Q, x . As an example 
of halo fluxes, Fig. 3 shows the solar-modulated flux of anti-protons from neu- 
tralino annihilations in the halo versus the neutralino relic density, together 
with the BESS 1-a measurement [12] , both at 1 GeV kinetic energy and solar 
minimum. The envelope of maximum fluxes decreases linearly, not quadrati- 
cally, with decreasing density. This shows that even if an increase in the cross 
section compensates the decrease in one of the powers of the density, fluxes still 
decrease linearly with the halo WIMP density. Hence WIMPs which are a sub- 
dominant component of the total CDM would be extremely difficult to detect 
through this method, as already concluded in [4, 5]. However, we expect that 
the intensity of the high-energy neutrino emission from the Sun and the Earth 
would in many cases remain high with decreasing density 

As the Sun and Earth orbit about the center of the galaxy, they sweep 
through the dark matter halo. Interactions with nuclei within the Sun and 
Earth slow WIMPs enough such that the WIMPs can become gravitationally 
captured. The WIMP capture rate of a macroscopic body is given by [13] 

c - ^ F (°) E ^Er 1 ^ «i) - <*.> p*> ( 4 ) 

where p x is the halo density of the WIMPs, m x is the mass of the particular 
WIMP, (a x _i) s is the elastic scattering cross section of the WIMP with the 
i-th type of nucleus in the body, /3j is the fraction of the body's mass in that 
nucleus, M is the total mass of the body, i s the squared escape velocity 

and (v^ ut 4 ) is the squared speed at which capture is kinematically cutoff, both 
averaged over the mass distribution of the i-th element. F(0) is the zero- velocity 
phase space distribution of the WIMPs and we call (a s ) the averaged scattering 
cross section appearing in the capture rate C. Furthermore, the annihilation 
rate Ta of captured WIMPs within the body can be written in terms of the 
capture rate as [14] 

F, = ^tanh 2 (£) (5) 

where t is the age of the body and ta — (CC'a)^ 1 ^ 2 is the equilibration time 
between capture and annihilation. Ca, the annihilation rate per WIMP pair, 
can be written in terms of effective volumes V\ and Vi as Ca = (<rv) A V2/V1, 
where the effective volumes depend on the pressure and temperature within 
the body in question (Vj = 3T n /(2jm x Gpo), where T and p are the core 
temperature and density of the Sun or Earth) . 

When the equilibration time is much shorter than the age of the body, ta <C 
t, the WIMP population is in equilibrium, i.e. the annihilation rate and the 
capture rate are equal (apart from a factor of two). In this case, 

C 

T A ^ — - {<Js)p x - (6) 
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If a decrease in WIMP relic density is compensated by a corresponding increase 
in the scattering cross section, the annihilation signals from Sun and Earth 
remain constant as the relic density decreases. 

If instead the equilibration time is much longer than the age of the body 
ta 3> t, then 



ri4 - C A Ch> ~ (a a v) (a s fpl ~ (a a v) J^-. (7) 



Thus if the scattering cross section rises as the annihilation cross section in- 
creases to decrease the current density of neutralinos, the annihilation rate 
may remain high. In this regime, overcompensation, namely a rise in the neu- 
trino flux despite decreasing relic neutralino densities, is possible. Namely, if 
(<7 S ) / {<r a v) E u is constant, then Ta ~ (c a v) ~ (a a v) E v ~ ft~ - This overcom- 
pensation is different from the increase in rates in direct and indirect detection 
qualitatively put forward in Ref. [5] (see Fig. 2 in the first paper of Ref. [5]), 
which for indirect detection would hold in the case of equilibrium between cap- 
ture and annihilation. Ref. [5] stated that rates either remain constant or in- 
crease slightly with decreasing densities. Below we see that this statement holds 
only at relatively high subdominant relic densities, while for smaller densities 
rates decrease (see Fig. 5 below for indirect detection, and Fig. 2 in Ref. [4] for 
direct detection). 



Notice that t a = {CC A ) 1/2 is 



{<? a v) 



-1/2 



TA ~ (V,) / \ ■ (8) 

Thus, if both (cr a v) and (a a v) E v increase about equally, then ta ~ (c s ) 1 ^ 2 . 
Hence a smaller ta corresponds to a larger scattering cross section. Since cap- 
ture is more efficient for large (a s ), the highest fluxes will come from the models 
with the smallest ta- 

Figs. 4a and 4b illustrate this clearly 1 . There we plot the muon rates 
in underwater or under-ice detectors (such as AMANDA [15], IceCube [16], 
ANTARES [17], and NESTOR [18]) from annihilation in the Sun and the Earth 
as functions of the equilibration time ta- (The vertical dashed line in figs. 3a 
and 3b denotes the age of the solar system, 1.5 x 10 17 s.) These muons are 
produced by the neutrinos generated in neutralino annihilations in the Earth 
or the Sun. These neutrinos interact with matter in or near the underwater or 
under-ice detectors through the charged current interaction vN — ► fiX produc- 
ing muons. In this paper, we take the energy threshold of detectable muons to 

1 Notice that while in the other figures we use a regular grid of points covering the region 
with models, in Figs. 4, 7, 8 and 9c, d and e below we show the original points in the table of 
models. No meaning should be assigned to the density of points, as the latter is an artifact 
of the scanning used to generate the table of models (see discussion in Ref. [6]). 
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be 25 GeV, which is a reasonable value for IceCube, and assume that it will be 
possible to detect a flux of 10 upward-going muons/km 2 /yr. 

For the Sun (Fig. 4a), ta can be <C t, ~ t or 3> t. The highest fluxes corre- 
spond to ta <C t, namely to models for which equilibration between capture and 
annihilation has taken place. So for the Sun, we do not expect overcompensa- 
tion. In fact, in Fig. 5, which shows the neutrino-induced up-going muon fluxes 
as a function of the relic density fl x h 2 , we see compensation down to densities of 
VL x h 2 ~ 10~ 2 , and a suppression of the rate for lower densities. In this figure we 
include all neutralino masses, since the flux of neutrinos from the Sun does not 
change much with neutralino mass. Present bounds from MACRO [19], Bak- 
san [20], and Super-Kamiokande [21] exclude fluxes higher than approximately 
10 3 muons/km 2 /yr. Models with Vt x h 2 ~ 1 x 10~ 3 may have signals above 10 
muons/km 2 /yr. So a kilometer-size detector may be able to probe even neu- 
tralinos which constitute only 0.7% of the dark halo using neutrino signals from 
the Sun. 

The Earth is either in the regime where ta 3> t or in the regime in which 
ta — t (see Fig. 4b). Thus the highest fluxes correspond to ta — t. Equilibra- 
tion between capture and annihilation is not reached, and so we expect cases 
with overcompensation. In fact, this is apparent in Fig. 6a, where we plot the 
muon fluxes from neutralino annihilations in the Earth as function of Q, x h 2 for 
neutralinos lighter than 150 GeV. The envelope of the highest fluxes increases 
with decreasing densities down to values of tt x h 2 ~ 3 x 10~ 3 . Close to this value 
of fl x h 2 , which corresponds to a neutralino fraction in the halo of only 2%, the 
muon fluxes become higher than 10 muons/km 2 /yr, while signals from the Earth 
from neutralinos constituting the whole of the halo would only reach the level 
of 1 event/km 2 /yr (with a threshold of 25 GeV). With the energy threshold of 
25 GeV, neutralinos lighter than 150 GeV are expected to be within reach of a 
kilomctcr-sizc neutrino telescope searching for neutrino signals from the Earth 
only if they are subdominant. 

For neutralinos heavier than 150 GeV, the envelope of the highest muon 
fluxes from the Earth remains almost constant with decreasing densities down to 
tt x h 2 ~ 5 x 10~ 3 (see Fig. 6b). There are expected rates above 10 muons/km 2 /yr 
for densities as low as Q x h 2 ~3x 10~ 3 , i.e. a kilometer-size neutrino telescope 
may detect neutrino signals from the Earth for neutralinos with a halo fraction 
from 1 to as low as 2%. 

The predicted fluxes from the Earth arc all below the current upper limits 
of approximately 10 3 muons/km 2 /yr from MACRO [19], Baksan [20], Super- 
Kamiokande [21] and AMANDA [22]. Notice that in similar figures appearing 
in the literature (e.g. in Ref. [22]) there arc models with signals higher than 
shown here. The reason is a lower experimental energy threshold (e.g. 1 GeV) 
and a different choice of the minimum density value the neutralinos can have to 
constitute the whole of the halo (Sl x h 2 = 0.025 instead of 0.15). 

If detection of a direct or indirect dark matter signal occurs, some model 
discrimination could be achieved by looking for a signal in other ways. This is 
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illustrated in Figs. 7a, 7b and 8a, 8b. 

In Fig. 7a we present the fluxes of up-going muons produced by neutrinos 
from ncutralino annihilations in the Earth versus relic density, for all models 
with fluxes from the Sun larger than 10 muons/km 2 /yr. This figure shows that 
only some of the models detectable from the Sun are also visible from the Earth. 
Contrariwise, all the models detectable from the Earth produce signals also 
detectable from the Sun, as can be seen in Fig. 7b, where we present the fluxes 
from the Sun of models with fluxes from the Earth above 10 muons/km 2 /yr, 
again versus relic density. Notice that the points plotted have densities ranging 
from 100% to 1% of the halo dark matter density. The big range of fluxes in 
Fig. 7a is due to the fact that the capture rate in the Sun depends on both 
the spin-independent and spin-dependent scattering cross sections, while the 
capture in the Earth only depends on the spin-independent scattering cross 
section. 

The points in Fig. 7b visible from both the Earth and the Sun in kilometer- 
size neutrino telescopes are also plotted in Fig. 8a, which shows that all of them 
are within the sensitivity reach of the current generation of direct dark matter 
detectors (such as CRESST II [23], ZEPLIN II [24], CDMS-Soudan [25], and 
EDELWEISS II [26]). Fig. 8a presents the spin-independent scattering cross sec- 
tion of neutralinos with protons multiplied by the neutralino halo fraction (the 
product that enters the direct detection rate) versus the neutralino mass, for 
models producing fluxes higher than 10 muons/km 2 /yr from both the Earth and 
(therefore also) the Sun. In the figure we include the CDMS [27] and EDEL- 
WEISS 2000-2002 [28] direct dark matter search limits, as these bounds are 
currently the most stringent. The CDMS and EDELWEISS 2000-2002 limits, 
which we have not imposed on the models, eliminate several neutralino models 
with high fluxes from both the Earth and the Sun. However, many models 
fall outside the current bounds but within the discovery limit of the current 
generation of detectors. 

Regarding indirect signals from the Sun, on the other hand, not even the 
larger direct-search detectors of the next generation (such as ZEPLIN IV [24] , 
Cryo Array [29], GENIUS [30], and XENON [31]) will be able to examine all 
models that give a detectable signal from the Sun in a kilometer-size neutrino 
telescope. In fig. 7b, the models producing more than 10 muons/km 2 /yr from 
the Sun are plotted in the rescaled spin-independent scattering cross section 
versus mass plane. Although several models fall within the sensitivity reach of 
future direct detectors, several others fall outside. The latter are those with 
small spin-independent cross sections but large spin-dependent cross sections 
with protons in the Sun. 

Independently of indirect signals, direct searches may reach much lower den- 
sities. Figs. 9a to 9e show all the models in our table separated in decades of 
relic density Q x h 2 or cquivalcntly halo density fraction fee dm- In Fig. 9a neu- 
tralinos constitute between 100% and 10% of the halo, in Fig. 9b between 10% 
and 1%, in Fig. 9c between 1% and 0.1%, in Fig. 9d between 10~ 3 and 10~ 4 and 
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in Fig. 9e between 1CP 4 and 1CP 5 of the halo (the smallest density fraction being 
2 x 10~ 5 ). In all these figures we plot the rescaled spin-independent scattering 
cross section (i.e. the spin-independent scattering cross section of neutralinos 
with protons multiplied by the neutralino halo fraction, which is the product 
that enters the direct detection rate) versus mass. We also include the present 
direct-detection bounds and future sensitivity limits. We clearly see in these 
figures that current data are probing neutralino halo fractions down to 1% (see 
Figs. 9a and 9b). In the next few years, data are expected from detectors such 
as CRESST II [23], ZEPLIN II [24], CDMS at the Soudan mine [25], and EDEL- 
WEISS II [28]. These experiments will probe halo fractions as small as 10~ 4 
(see Figs. 9a,b,c,d). Further in the future, another generation of detectors, such 
as ZEPLIN IV [24], CryoArray [29], GENIUS [30], and XENON [31], may be 
able to reach down to neutralino halo fractions of order 10~ 5 , which are the 
smallest in our table of models. 

In summary, indirect detection rates of neutralinos annihilating in the Sun 
and the Earth remain above present detection thresholds even for neutralinos 
constituting only 1% of the halo dark matter. Indirect rates are drastically re- 
duced at lower density fractions. We find one instance in which the detectability 
is highest around a halo fraction of 1%, namely the case of light neutralinos giv- 
ing neutrino signals from the Earth. We do not find that the detectability of relic 
neutralinos using signals from Sun or Earth is usually favored for neutralinos of 
small Q x h 2 , as was previously stated [5]. 

Present limits from direct detection experiments probe down to 1% of the 
halo density, while data expected in the near future will probe down to 10~ 4 
and the next generation of detectors may reach 10~ 5 of the halo density (as 
already discussed in Ref. [4]). 

Many neutralino models constituting from 1% to 100% of the dark matter 
halo are detectable in three ways (direct, indirect from the Sun, indirect from 
the Earth), some in two, and some in only one. Thus if a signal is found either in 
direct detection experiments or in indirect searches for neutralino annihilations 
in the Sun or Earth or in both, the question of which component of dark matter 
was found, the primary or a sub-dominant one, may remain open. However, if a 
signal would be seen from annihilation in the halo or in the center of the galaxy, 
it would strongly point towards neutralinos being the dominant component of 
the dark halo, since the signal from annihilation in the halo decreases sharply 
with halo fraction. 
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Figure Captions 



Fig. 1 Spin-indcpcndcnt ncutralino-proton scattering cross section cr x _ p , zero 
momentum (and zero temperature) limit of the neutralino-neutralino anni- 
hilation cross section times relative velocity {(J a v), neutralino relic density 
il x h 2 , and Q x h 2 if coannihilations would not occur ("w/o coanns"), as 
functions of the gaugino fraction Z g of the neutralino for a specific "line" 
of models with fixed neutralino mass and all other MSSM parameters held 
constant (see text). 

Fig. 2 Neutralino relic density Q, x h 2 versus the inverse of the zero momen- 
tum (and zero temperature) limit of the neutralino-neutralino annihilation 
cross section times relative velocity {a a v) 1 together with the approximate 
relation in Eq. (2) (straight line) . A regular grid of points shows the region 
covered with models. 

Fig. 3 Solar modulated anti-proton flux at 1 GeV kinetic energy at solar mini- 
mum (using the force-field approximation as implemented in DarkSUSY) 
as a function of the neutralino relic density. A regular grid of points shows 
the region covered with models. The horizontal band shows the BESS 1-cr 
measurement of the antiproton flux. 

Fig. 4 Flux of muons in underground detectors with > 25 GeV from neu- 
trinos produced in neutralino annihilations in: (a) the Sun, (b) the Earth, 
as a function of the equilibration time in the Sun or Earth, respectively. 
Each point represents an actual model. The dashed line shows the current 
age of the Solar system, 1.5 x 10 17 seconds. 

Fig. 5 Flux of muons in underground detectors with > 25 from neutrinos 
produced in neutralino annihilations in the Sun as a function of the neu- 
tralino relic density. A regular grid of points shows the region covered 
with models. Approximate indirect detection current bounds and future 
sensitivity of a kilometer-size detector are indicated. 

Fig. 6 As Fig. 5, but for neutralino annihilations in the Earth and neutralinos 
with: (a) masses m x < 150 GeV, (b) masses m x > 150 GeV. 

Fig. 7 Flux of muons in underground detectors from neutrinos from neutralino 
annihilations in: (a) the Earth, (b) the Sun, as a function of the neu- 
tralino relic density, for models which concurrently have muon fluxes from 
neutralino annihilations in (a) the Sun or (b) the Earth with rates greater 
than 10 muons/km 2 /yr. Each point represents an actual model. 

Fig. 8 Spin-independent neutralino-proton cross section times neutralino halo 
fraction as a function of neutralino mass, for models which (a) have muon 
fluxes from neutralino annihilations both in the Earth and in the Sun with 
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rates greater than 10 muons/km 2 /yr, and (b) have muon fluxes from the 
Sun with rates greater than 10 muons/km 2 /yr. Each point represents an 
actual model. Best direct detection current bounds and selected future 
sensitivity limits are shown. 

Fig. 9 Spin-indcpcndcnt scattering cross section of neutralinos with protons 
times neutralino halo fraction versus mass, together with present direct 
detection bounds and future sensitivity limits. Models are separated in 
decades of halo fraction fee dm '■ (a) 0.1 < Jccdm < 1; (b) 10~ 2 < 
fccDM < 0.1; (c) 10- 3 < fccDM < 10- 2 ; (d) 10- 4 < fccDM < HT 3 ; (e) 
10~ 5 < fccDM < 10~ 4 (the smallest value we found is Jccdm — 2 10 -5 ). 
In Figs. 9a and 9b, a regular grid of points show the region covered with 
models, in the others each point represent an actual model. 
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